Average entanglement for Markovian quantum trajectories 
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We study the evolution of the entanglement of noninteracting qubits coupled to reservoirs under 
monitoring of the reservoirs by means of continuous measurements. We calculate the average of 
the concurrence of the qubits wavefunction over all quantum trajectories. For two qubits coupled 
to independent baths subjected to local measurements, this average decays exponentially with a 
rate depending on the measurement scheme only. This contrasts with the known disappearance of 
entanglement after a finite time for the density matrix in the absence of measurements. For two 
qubits coupled to a common bath, the mean concurrence can vanish at discrete times. Our analysis 
applies to arbitrary quantum jump or quantum state diffusion dynamics in the Markov limit. We 
discuss the best measurement schemes to protect entanglement in specific examples. 

PACS numbers: 03.67.Pp, 03.65.Yz, 03.67.Mn 



I. INTRODUCTION 



Entanglement is a key resource in quantum informa- 
tion. It can be destroyed or sometimes created by inter- 
actions with a reservoir. When the two non-interacting 
parts of a bipartite system are coupled to independent 
baths, entanglement typically disappears after a finite 
time [iHJ|. This phenomenon, called "entanglement sud- 
den death" (ESD), occurs for certain initial states only or 
for all entangled initial states, depending on whether the 
system relaxes to a steady state belonging to the bound- 
ary of the set of separable states (e.g., to a separable pure 
state for baths at zero temperature) or to its interior (e.g., 
to a Gibbs state at positive temperature) @ . A quantum 
state Ues on this boundary if it is separable and an arbi- 
trarily small perturbation makes it entangled; this is the 
case, for example, for a pure separable state. When the 
two parts of the system are coupled to a common bath, 
sudden revivals of entanglem_ent_may take place after the 
state has become separable 



In this article we consider the loss of entanglement be- 
tween two non-interacting qubits coupled to one or two 
baths monitored by continuous measurements. Because 
of these measurements, the qubits remain at all times in a 
pure state \tjj{t)), which evolves randomly. To each mea- 
surement result (or "realization") corresponds a quan- 
tum trajectory t G K-|. i-)- |V'(i)) in the Hilbert space 
of the qubits. In the Born-Markov regime, the dynamics 
is given by the quantum jump (QJ) model d, (To] or, in 
the case of honiodync and heterodyne detections, by the 
so-called quantum state diffusion (QSD) models [10l - [l2| . 
We study how the entanglement of the qubits evolves 
in time by calculating the average C^(t) of the Wootters 
concurrence of \ip{t)) over all quantum trajectories; C^(t) 
differs in general from the concurrence Cp(t) of the den- 



sity matrix p{t) = \^it)){'4>(t)\ (here and in what follows 
the overline denotes the mean over all quantum trajec- 
tories) [isl. For two qubits coupled to independent 
baths, we find that 



(1) 
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where Co = C^(o) is the initial concurrence and k > de- 
pends on the measurement scheme but not on the initial 
state In particular, if Cq > and iESD S]0, oo[ is 

the time at which entanglement disappears in the density 
matrix (assuming that this time is finite), then Cp(t) = 
at times t > tssD whereas C^(t) can only vanish asymp- 
totically. The continuous measurements on the two baths 
thus protect on average the qubits from ESD. Of course, 
this does not mean that all random wavefunctions |'0(i)) 
remain entangled at all times. But in some cases, such as 
for pure dephasing or for infinite temperature baths, one 
can find measurement schemes such that k = 0; then, 
for all trajectories, if the qubits are maximally entangled 
at t = they remain maximally entangled at all times. 
We show that the best measurement scheme to protect 
entanglement is in general given by homodyne detection 
with appropriately chosen laser phases. Related strate- 
gies using quantum Zcno effect 11 51. entanglement dis- 
tillation Il6[ , quantum feedback [17| . and encoding in 
qutrits [l8| have been proposed. It is assumed in this 
work that the measurements on the baths are performed 
by perfect detectors. The impact of detection errors has 
been studied in [isj . 

When the qubits are coupled to a common bath, we 
find that C^(t) has a more complex time behavior than in 
([1]). It may vanish at finite discrete times, and, for some 
initial states, be equal to Cp(t) . It is worthwhile to stress 
that the formula ([1]) is valid provided not only each qubit 
is coupled to its own bath, but also the baths are mon- 
itored independently from each other by the measure- 
ments. This means that the measurements are performed 
locally on each bath. Instead of looking for the measure- 
ment scheme maximizing the average concurrence C^(t) 
of the two qubits in order to obtain the best entanglement 
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protection, it is also of interest to find a way to perform 
the measurements such that C^(t) is minimum and co- 
incides with the concurrence Cp(i) of the density matrix. 
This problem has been studied numerically in Ref. 
and analytically in [l9| for specific models of couplings 
with the two baths. Our result ([T]) implies that for any 
Markovian dynamics, if the two qubits are initially en- 
tangled and ESD occurs for the density matrix p{t), a 
scheme with the aforementioned property must necessar- 
ily involve measurements of non-local (joint) observables 
of the two baths. In the models studied in non- 
local measurements are indeed used in order to obtain an 
optimal scheme satisfying C^(t) = Cp(t) ■ 

The paper is organized as follows. We briefly recall 
in Sec. the definition of the concurrence of pure and 
mixed states and review the quantum jump unraveling 
of a Lindblad equation for the density matrix in Sec. IIIII 
We treat the simple and illustrative case of two two-level 
atoms coupled to independent baths at zero temperature 
in Sec. |lVl before showing formula ([1]) in Sec. |V] for a 
general quantum jump dynamics. The QSD unravelings 
are considered in Sec. IVIl we obtain the average concur- 
rence for such unravelings as limits of the concurrence 
for QJ dynamics (corresponding to homodyne and het- 
erodyne detections with intense laser fields) . Section IVIII 
is devoted to the evolution of the entanglement of two 
qubits coupled to a common bath at zero temperature. 
The main conclusions of the work are given in Sec. IVIIII 

II. ENTANGLEMENT MEASURES FOR 
QUANTUM TRAJECTORIES 

The entanglement of formation of a bipartite quan- 
tum system 5 in a pure state \ip) is defined by means 
of the von Neumann entropy = — tr(/9A InpA) = 
— tT{pB^^ Pb) of the density matrices pA ~ tTB{\ip){'4'\) 
and pb = tr^ (I?/') (t/jI) of the two subsystems A and 
B composing S [20|. If 5 is in a mixed state, Ep is 
the infimum of '^PkE^Pf, over all convex decompositions 
p = J2k Pkl'ipk) (^fc I of its density matrix (with pk > and 
\\tpk\\ = !)• When A and B have two-dimensional Hilbert 
spaces, Ep = f{Cp) is related to the concurrence [2l| Cp 
by a convex increasing function / : [0, 1] — >■ [0, ln(2)]; p is 
separable if and only if Cp ~ 0, i.e., Ep = 0. For a pure 
state 

C4, = \{ay®ayT)^,\ (2) 

where ay = i(|i)(t| - It) (ID is the y-Pauh matrix, T : 
IV") = J2s,s' ^ss'\s,s') Y.^ ,,, c*^,|s, s') the anti-unitary 
operator of complex conjugation in the canonical basis 
{\s,s') = \s) \s');s,s' =t,i} of C2 ® C2, and (•)v, = 
■ \ip) the quantum expectation in state {-ip). 
For quantum trajectories, one has always i?(/>(t) > 
Ep(^t) , this inequality being strict excepted if the decom- 
position 



realizes the infimum defining Ep(^ty Thanks to the con- 
vexity of /, > /(C'v'(*))' Thus equation (jT]) shows 
that for independent baths and if Co > 0, £'^,(t) > 
/(Coe^'**) > whatever the measurement scheme. 

It is legitimate to ask which entanglement measure 
should be averaged since, for example, £'^,(t) = Eq could 
be constant and C^(t) time-decreasing if Eq ^ 0,ln2. 
The concurrence is a natural candidate as it corresponds 
for pure states to the supremum over all self-adjoint local 
observables Ja and Jb with norms less than one of the 
modulus of the correlation between J a and Jb , 

C^(t) = sup \{Ja® JB)4,(t} 

\\JaI\\Jb\\<1 

-(Ja ® lB)^0(t)(lA «i Js)0(t)| • (4) 

Moreover, C^,(t) is easy to calculate in the Markov regime 
and gives a lower bound on E^(ty 

III. QUANTUM JUMP MODEL 

Let us briefly recall the Q J dynamics [1, [12, . As a 
result of a measurement on a particle (e.g. a photon) of 
the bath scattered by the qubits, the qubits wavefunction 
suffers a quantum jump 

where the jump operator is related to the particle- 
qubits coupling and the indices m, i label all possible 
measurement results save for the most likely one, which 
we call a "no detection". In the weak coupling limit, the 
probability that a measurement in the small time interval 
[t^t + dt] gives the result (m, i) is very small and equal 
to dpj,,(i) = 7^11 J^|V'(i))pdt. The jump rate does 
not depend on \ip{t)) and is proportional to the square of 
the particle-qubit coupling constant. In the no-detection 
case the wavefunction of the qubits evolves as 

+ - ||e-^^e..^*S^))r ''°" = ''° 2 ^ ^-^^-^^ 

(6) 

where Hq is the Hamiltonian of the qubits. The proba- 
bility that no jump occurs in the time interval [tQ,t] is 
Pnjito,t) = ||e-'^-«(*^*°)|V'(to))|P- (This is proven as fol- 
lows: a.S Pnj{tQ,t)-pnj{to,t + dt) = Em,idp^(t)Pnj(to,i), 

one has In p^^ (to, t)/dt = - E™. tLII^^IV'W)!!' = 
(9/9t)ln|le-'^<=«(*-*«)|i/'(io))|P by ©.) It is not difficult 
to show that the density matrix p{t) = |?/'(t)) (V'(OI 
satisfies the Lindblad equation 

^ = -i[Ho,p]+J2 f^[ji,,pJii - \{JilJLp}) (7) 

m,2 



p{t) = m)){m\^ I dp[i:mt)){m\ (3) 



where {•,•} denotes the anti-commutator. It is known 
that many distinct QJ dynamics unravel the same master 
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equation ([7]) [22] ■ For two qubits coupled to independent 
reservoirs Ra and Rb, the jump operators are local, i.e., 
they have the form 

J^^Ib , lA®J,n (8) 

depending on whether the measurements arc performed 
on Ra or Rb- Here are 2x2 matrices. 

The aforementioned absence of ESD for the mean con- 
currence of two qubits coupled to independent baths can 
be traced back to the existence of trajectories for which 
\ilj{t)) remains entangled at all times. Actually, for a tra- 
jectory without jump, |V'nj(i)) « e^'*-^"'* |V'(0)), see ([6]). 
By (jHl) and since the qubits do not interact with each 
other, 6^'*^°''^ is the tensor product of two local opera- 
tors acting on each qubit. If |V'nj(0) would be separable 
at a given time t then, by reversing the dynamics (i.e., by 
applying e'*^"'' to |V'nj(0)) ^'^'3 would deduce that \ip{0)) 
is separable. Hence C'^„j(t) > if Co > 0. But the no- 
detection probability between times and t is nonzero 
and thus C^(t) > at all times. Note that this argument 
does not apply if non-local observables of the two baths 
are measured or if the two qubits are coupled to a com- 
mon bath, since then the jump operators are non-local. 

IV. PHOTON COUNTING 

Let us illustrate the random dynamics described pre- 
viously on a simple and experimentally relevant exam- 
ple [2^. Each qubit is a two- level atom coupled reso- 
nantly to the electromagnetic field initially in the vac- 
uum (zero-temperature photon bath). The two atoms 
are far from each other and thus interact with indepen- 
dent field modes. Two perfect photon counters Di make 
a click when a photon is emitted by qubit i {i = A, B) , 
whatever the direction of the emitted photon. Doing 
the rotating wave approximation, the jump operators are 
Ji = ai_ = |i)(tl- For simplicity we take Hq = 0. 
By ([6]), if no photon is detected in the time interval [0, t] 
the qubits state is 

1^(0) =AA(i)-i ^ c,,,e-^»»'*/2|5,s') (9) 

with = lA+lB, It I = 7A, lit = 7B, 7ii = (7^ 
being the jump rate for detector Di), Css' = {s, s'\Tp{0)) , 
and N{t)^ = Icss'pe"'^-'*. The concurrence © 

of |?A(0) is C{t) ^ CoAA(t)-2e-(TA+7B)t/2 ^jth Cq = 
2|c-t-|-C4,4, — c-t'4,C4,-|-|. If a photon is detected at time tj 
by, say, the photon counter Da, the qubits are just 
after the jump ([5]) in the separable state |V'(^j+)) ^ 
\i) (g> (ctte-''ttt372||) + ct4.e-'^t.tt372||)). Since neither 
a jump nor the inter-jump dynamics can create entangle- 
ment (the jump operators ^ being local), |'0(t)) remains 
separable at all times t > tj, even if more photons are 
subsequently detected. Thus C{t) = if at least one pho- 
ton is detected in the time interval [0,t\. Averaging over 
all realizations of the quantum trajectories and using the 
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FIG. 1: (Color online) Concurrences of two qubits coupled to 
independent baths at positive temperature as a function of 
for 7; = 7i/2 = 7 and |V(0)) = (|tt) - m)/V2: (2a) C.^t) 
for the density matrix (blue dashed line); (2b) C^{t) for a 
single trajectory (black dotted line); (2c) C^(t) averaged over 
1500 trajectories and from Eq. ([1} (red solid lines); (2d) C^(t) 
for the best measurement scheme (see text). 

probability pnj(0,t) = Af{t)'^ that no photon is detected 
in [0,t], one finds C{t) = Coe-(T-*+''«)*/2^ 

This argument is easily extended to baths at positive 
temperatures by adding two jump operators J\_ = a\_ 
with rates < 7L. The mean concurrence is then 

Cjt) = C'oe-('^^+'^-+''++''^'*/^ It is compared in Fig. 
□ with the concurrence of the density matrix obtained by 
solving the master equation ([7]), which shows ESD for all 
initial states. 



V. GENERAL QUANTUM JUMP DYNAMICS 

We now consider a general QJ dynamics with jump 
operators given by ([8]). The Hamiltonian of the qubits 
has the form Hq = Ha (8) 1_b + 1a ^ Hb- Let K = 

KA<E)lB + iA<i^KB with 

K^ = lY.^rnJ:U^, (10) 

rn 

7*„ being the jump rates for the detector Dj [i = A,B). 
We first assume that no jump occurs between t and t+dt. 
By expanding the exponential in one gets 

C{t + dt) = p„j(t, t + di)-i I {ay ® OyT)^^^^ (11) 

+-iAt{Hl^ay (g) UyT ^Oy® ayTH^s)^^^.^ + 0{dtf\ 

where pnj(t,i + d<) = {l-2Kdt + 0{dty)^(^t) is the prob- 
ability that no jump occurs between t and t + dt. Now, 
for any local operator Oi acting on qubit i, one has 

{0^<7y®<7yT)^^^^^{<7y®<JyT0l)^^^ =^tTc2m (l2) 
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with 



C{t) = {ay ® ayr).^(t) =2(4^(t)cJ^(t) - 4^(t)cJ^(t)) 

(13) 

and Css'{t) = {s,s'\ilj{t)). Since C(i) = |C(i)| and iJeff = 
Y,t{Hz - iA'j), this gives 



C(t + dt)p„j(i,t + dt) = C(f)(l-trc4(if)^ + 0(dt2)) . 

(14) 

If detector Di gives the resuh m in the time interval 
[i, i + di], the concurrence is by virtue of (O 



d™^;^ + dt) = 3%^C(i) I detc. iJi^)\ 



'jump ' 



where we have used the identity 

{Ojay ® <JyTO^)^,(t) - C(t)detc2(0j) 



(15) 



(16) 



vahd for any local operator Oi acting on qubit i. Collect- 
ing the previous formulas and usi ng the M ar kov p rop- 
erty of the jump process, one gets C{t + dt) = C{t){l — 
KQj dt -^^(dt^)) with 



KQJ = ^trc4(K) -^7;,|detc2(j; 



(17) 



Letting d< go to zero, one obtains dC(t) /dt = — kqj C{t). 
The solution ([1]) of this differential equation has the 
exponential decay claimed previously. To show that 
'tQj ^ 0; let 20^ be the argument of dctca ( J^). We write 

«QJ = l:ra,^lln[^^C-{JM - 25R{e-2i«™ dctc^ ( J^,)}] /2 

as 



E¥(Kti4Jt)-(;iJ^i;)| 



|(t|25R4j;)| 



(18) 

with e-'**" and 23?J;„ = 4, + J^t. Thus kqj is 

non-negative. 

Note that kqj = if all matrices are self-adjoint 
and traceless (then 9)^ = 7r/2 and U.P^ = 0). We 
show in Fig. [2] the concurrence of the density matrix 
given by solving ^ for a pure dephasing with J' ~ 
gi7r/4^i_ _^ g-i7r/4^j^^ Qnc has BSD for all initial states 

save for^O)) = (|n) ± Since kq,, = 0, © 

implies C(t) = Cq. If the two qubits are maximally en- 
tangled at t = 0, then C^(t) = C(t) = Co = 1 for a// 
quantum trajectories at any time t > 0. Therefore, for 
pure dephasing one can protect perfectly the qubits by 
measuring continuously and locally the two independent 
baths. 

We can now give the optimal measurement scheme to 
protect the entanglement of two qubits coupled to in- 
dependent baths at positive temperatures. Let us re- 
place the photon-counting jump operators Jj_ = tr!^ by 
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FIG. 2: (Color online) Same as in Fig[T] for pure dephasing 



and the initial state [^(O)) = ^(|tt) + e~"^ii|)): (la) Cp^t) 
for if 



^ (blue dashed line); (la') Cp{t) for = Q (blue line 
showing ESD); (lb,lc) C^(t) — C^i^t) (red solid line). 



are unitary 2 x TV matrices. This corresponds to a rota- 
tion of the measurement basis and gives another unravel- 
ing of the master equation ((T]). Let us stress that the new 
jump operators J* still act locally on each qubit. By (fTT)) . 

the new rate is k = Eu,^(\n^l<-|- \ml<+l)V2- By 



using = 1 and optimizing over all unitaries Ui, 

one finds that the smallest disentanglement rate arises 
when, for example, = iu'j^ = l/\/2 (N = 2) and is 
given by 



■'qj 



(19) 



i=A,B 



Note that Kq^* 



kqj at zero temperature and Kqj' ~ 
(perfect protection) at infinite temperature. The decay 
of C (t) for this optimal measurement is shown in Fig. [1] 
(green dashed-dottcd line). 



VI. HOMODYNE AND HETERODYNE 
DETECTION 

Let us come back to our example of two atoms cou- 
pled to the electromagnetic field initially in the vacuum. 
If homodyne photo-detection is used instead of photon 
counting, the jump operators become Jj_q, = cr!_ ± cti , 
being the amplitude of a classical laser field (there are 
now four jump operators since each homodyne detector 
involves two photon counters) [ll|. Assuming that the 
two photon beams emitted by the atoms are combined 
with the two laser fields via 50% beam splitters, the jump 
rates associated with Jj.^ are equal, 7j_^ = ji/2. Thanks 
to (fT7|) . one easily finds that the disentanglement rate for 
the new QJ dynamics, KQj(a) = {'Ja+Jb)/'^, is the same 
as for photon counting. 
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In contrast, KQj(a) depends on the laser amplitudes 
for pure dephasing (jump operators Jj_^ = • cr ± a; 
with Vj e M^, ||vi|| = 1, and a the vector formed by 
the Pauli matrices ax, <Jy, and ct^): then KQj(a) = 
2 7i min{a|, 1} for real a^'s. One reaches perfect 
entanglement protection {C{t) = Co) only for vanish- 
ing laser intensities of. In the case of two qubits cou- 
pled to two baths at positive temperatures, a general 
choice of jump operators such that the density matrix 
^ satisfies the master equation ([7]) with the four Lind- 
blad operators ctj., i = A,B, is ±q = ± 
with the jump rates 7* = lli/'^i laser amplitudes 



and Jl = T, 



bitrary unitary matrix (u^^)™^]^.. jy (see the discussion 
in the preceding section). The corresponding disentan- 
glement rate, KQj(a) = E^,, T^ftrc^ ( J// J/J + 2|a^|2 - 
2| detc2( J^) + (a;,)2|]/2, is equal to kqj(O) if det(J^) = 
or for complex laser amplitudes = ja^je'^'' satisfy- 
ing 20^ = arg(det(J^)); otherwise, AtQj(a) is larger then 
'^Qj(O)- We can conclude that the smallest disentangle- 
ment rate is given by (jl9p and the best unravelings to 
protect the entanglement of the qubits are either the Q J 
model with jump operators J{ oc (7!).) 2CT^-t-(7^) scr^ and 

J2 cx (7^)5ct!^ — (7!_)^cr!_ or the corresponding homodyne 
unraveling with laser phases 9\ = 7t/2 and 62 = 0. 

Let us now consider a general QJ model with jump 
operators J^. A new unraveling of (O is obtained from 
the QJ model with jump operators Jm,±a = ^m^'^m and 
rates 7m.±Q = Iml'^- For large positive laser amplitudes 
^ 1, this dynamics converges after an appropriate 
coarse graining in time to the QSD model described by 
the stochastic Schrodinger equation [111, 



|d^) = (-iifo - im + ^ /^(j;, - ^{.rj^ 



xd< + 7;„(3?(j; 



V.) (20) 



where dw^ are the Ito differentials for independent real 
Wiener processes satisfying the Ito rules dw^dw^ = 
5ij5mn'^i- One can determine the mean concurrence for 
the QSD model ([^0]) by taking the limit of the mean 
concurrence for the QJ dynamics with jump operators 



This gives again the exponential decay ([T]) but 



with a new rate 

«ho = -E7m(5?detc.(j;,) + l{'itTc2irjy 

. (21) 

In fact, if 26'J„ is the argument of det( J^^ ± aj,^) = 
(aj„)2 ± aj„tr(j;„) + 0(1) then for aj„ » 1, < > 0, 
one has e^'^^.±° ~ l±i5tr(j;,)/aj„. Using (Ull), a short 
calculation gives (pij) . 

Unlike kqj, Kho changes when the operators J^j in 
(PU]) acquire a phase factor, — >■ e"'^™ J^. This arises 
for homodyne detection with complex laser amplitudes 



'^m ~ l"^ml^'^"i I'^ml 1- Minimizing over the laser 
phases 9^^ yields 



opt 
«ho = 



i trc4 {!<) - 5^ 7™ (|det( J^) - 1 (trc. (4,))' 



2 



(22) 



It is easy to show that < kqj, this inequality be- 
ing strict excepted if the two eigenvalues of J,', have the 
same modulus for all (to,i). Thus optimal homodyne 
detection protects entanglement better than - or, if the 
aforementioned condition is fulfilled, as well as - photon 
counting. Let us stress that the optimal measurements 
(in particular, the laser phases 0^ minimizing the rate 
Kho) only depend on the Lindblad operators in the 
master equation ([7]) and are thus the same for all initial 
states of the qubits. 

Let us now discuss the case of heterodyne detection. 
The corresponding jump operators Jm.±a{^q) = ^ 
Q;J„e'^"*« depend on the time tq of the g-th jump due 
to the oscillations of the laser amplitudes [l^l- The as- 
sociated rates are 7m_±a = 7m/2 as for homodyne de- 
tection. We assume here that > 0. In the limit 
(a^)^ fll^/^l^ » 1 of large laser intensities and 
rapidly oscillating laser amplitudes, the QJ dynamics 
with jump operators Jm.±ai'^q) converges to the QSD 
model given by the stochastic Schrodinger equation [22 1 



-iHo - K)dt {^■^'rn)*^ ^™ 

|(j^)^|')dt + ^ V7^((j;,-i(j;,)^,)dc 

(^™>;(dC)*)llV'> (23) 



1, 
~V 

1 



where d^^ are the Ito differential of independent complex 
Wiener processes satisfying the Ito rules d^^^d^^j = and 
dCm(dC^i)* = SijSmndt. Eq. (f23l) describes the coarse- 
grained evolution of the normalized wavefunction |^/'(i)) 
on a time scale At such that (i) many jumps and many 
laser amplitude oscillations occur in a time interval of 
length At and (ii) \ip{t)) does not change significantly on 
such a time interval. These conditions are satisfied when 
("m)^7mAt > ni,^At > 1 and 7;„At < 1. We now 
show that the mean concurrence for the QSD model ([23]) 
is given by ([T]) and determine the rate k of its exponen- 
t ial de cay. This can be done by calculating the derivative 
dC{t)/dt in a similar way as in SeclYl using and the 
Ito rules. It turns out to be simpler to estimate directly 
the average concurrence of the QJ model for heterodyne 
detection in the aforementioned limits, in analogy with 
our previous analysis for homodyne detection. Let us 
first remark that the results of Sec. IVl remain valid if the 
jump operators Jm{t) vary slowly in time, on a time scale 
much larger than the mean time {am)~'^ /im be- 
tween consecutive jumps. Hence dC/dt = — Khct(0 C'(i) 
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and thus C{t) = Cq e~ ^' "^hotC*') ^ time-dependent 
rate Khet(i) given by ([18]). To simplify notations, we tem- 
porarily omit the sum in (fT5)l and do not write explic- 
itly the lower and upper indices m and i. Let us set 
T = tr(J)/2 = \T\e''^ and S = det(J) = e^'^\5\. Let 
26'±Q(t) denote the argument of det(J± ae'^*). Gener- 
alizing the calculation outlined above for homodyne de- 
tection, one gets e^f'i''^*) ~ e2'f"(l ± 2i3{r e~'"*}/a) as 
a > 1. By (Ull) this yields 



-2ie'^"3{Te-'"*}|^ 



J%)\ 



trc4(A') 



- 7(|(5| cos(26i - 2nt) + 2|rp sin2((p - nt)) 



up to terms of order a ^. 
latory integral /p dt' cos(20 



By neglecting the oscil- 
2r2i') (which is of order 
fl^^ <^ At < t) and approximating Jq dt' sin^ {tp — fit') 

by t/2, one obtains /J dt' KhGt(i') ~ t(trc4(X)/2-7|Tp). 
Putt ing together the previous results, this shows that 
C{t) -J> Coe""''"'* in the limit > $1/7 > 1 and 
n > (At)-i > 7, with 



^het 



trc4(i^) 



7m|tr(j;,)| 



(24) 



We note that ^hct ^ 



opt 



For given jump operators 



J^, the measurement scheme which better protects the 
qubits against disentanglement is thus given by homo- 
dyne detections with optimally chosen laser phases. In 
this scheme, the average concurrence decays exponen- 
tially with the rate ((22|) . 

Although ([23]) is different from the QSD equation for 
the normalized wavefunction introduced by Gisin and 
Percival [l^, the quantum trajectories t 1— \ip{t)) for 
the two dynamics are the same up to a random fluctu- 
ating phase [l^ which does not affect the concurrence 
C^(t). More generally, one can show that the mean 
concurrence for the QSD model with correlated com- 
plex noises satisfying the Ito rules d^^d^^ = w^„dt and 
d^m(d^ri)* — SijSmndt [5^, which givcs back the model 
of Gisin and Percival when uj^^ = 0, decays exponen- 
tially as in ([1]) if the two baths are independent, i.e., if 



,AB 



for any to, n. 



VII. QUBITS COUPLED TO A COMMON BATH 

We focus here on a specific model of two qubits 
with equal frequencies coupled resonantly to the same 
modes of the electromagnetic field initially in the vac- 
uum. A photon counter D makes a click when a pho- 
ton is emitted by qubit A or B. The jump opera- 
tor in the rotating wave approximation, J = ® 
1b + 1a ® cr_, is now non-local. We take — 












■1 t3c\ \ 




I \ 

> \ 









0,0 0.5 



Dimensionless Time yt 



FIG. 3: (Color online) Concurrence of two qubits coupled 
to a common bath versus yt for \ip{0)) = "^Iti) + "Tf l-l-t)- 
(la) Cp(t) (blue dashed line); (lb) C^(t) for a single trajec- 
tory (black dotted line); (Ic) C^(f) given by (|25p (red line 
superimposed on the blue line). Inset (2) is the same for 

i^(o)) = ^in) + ^ii;>. 



0. Proceeding as for independent baths, the contri- 
bution to the mean concurrence of quantum trajecto- 
ries without jump between and t is pnj(0, t)Cnj(t) = 
|((Ty (E) (Tyr)e-tff |^(Q)^ I and can be determined with the 
help of ([T3|) . By calculating the exponential oi K = 
7jtj/2, one finds e-(*-*°)^|V(to)) = E.,.' c.s' (ON, s') 
with ctt(t) = e-T(*-*«)ctt(to), 2css'{t) = (e^^(*-*o) + 
l)css'{to) + (e-T(*-*°) - l)cs'sita) for ss' =ti or if, 
and cii{t) = cii{tQ). Quantum trajectories having one 
jump in [0, t] give a nonzero contribution. The proba- 
bility density that the jump occurs at time tj € [0,t] is 
given by 7p„j(t„t)|| J|V'(t,-))|pp„j(0,t,) = 7 Mj,*, (t)' 
with 7Vij,t,(t) = ||e-(*-*^)^Je-*^-^|^/>(0))|| (this follows 

from the formula p„j (io, i) = ||e-(*-*°)^|V'(<o))P, see 
Sec. IIII[) . The contribution of trajectories having one 
jump in [0,t] is then obtained by multiplying this den- 
sity by Cij^tj (t) = SA/'ij,*^ (t)-2e-27t|c^^|2 ^nd integrating 
over tj. After two clicks, I'lpit)) = is in an invariant 
separable state. Therefore, trajectories with more than 
one jump do not contribute to the mean concurrence. 
Setting c± = c-f^ ± c^-f, one gets 



(25) 

The time behavior of the concurrence ([25]) depends 
strongly on the initial state. Unlike in the case of in- 
dependent baths, C{t) may vanish at nonzero finite dis- 
crete times to- A necessary and sufficient condition for 
this loss of entanglement (immediately followed by a re- 
vival) is c-i-t = and arg(ci-4.) = arg(c4,t-) (i.e., e+/ c- G 
] — oo, — 1[U]1, oo[). If this condition is fulfilled, C{t) van- 
ishes at time to = 7"^ ln(|c+/c_|), see Fig. [31 It is not 
difficult to show by solving the master equation ([T]) with 
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J = cr_ (g) Is + <8) cr_ that, for any initial state con- 
taining at most one excitation (i.e., such that c-t~|- = 0), 
C{t) = — c^e^^'''*|/2 coincides at all times with the 
concurrence Cp(t) for the density matrix. In contrast, if 

c^l 7^ then C{t) increases at small times whereas C'p(() 
decreases, as sh own in the inset of Fig. [31 For any ini- 
tial state, C{t) converges at large times i ^ to the 
same asymptotic value Coo = |c-P/2 as the concurrence 

Cpit) mm- 

A non-local measurement scheme depending on the ini- 
tial state \ip{0)) and such that C{t) = Cp[t) at a-H times 
t G [0,tEDs] has been found recently [l^ for two qubits 
coupled to two baths at zero temperature in the rotating- 
wave approximation. If one neglects the Hamiltonian of 
the qubits, this scheme is time-independent. The corre- 
sponding quantum trajectories are given by a QSD equa- 
tion (26| for homodync detection with two jump operators 
Ji and J2 similar to the jump operator J introduced in 
this section, combined with intense laser fields via 50% 
beam splitters, as described in Sec. |VI] (the main dif- 
ference between J1.2 and J comes from the presence of 
appropriately chosen phase factors in front of a_ and a+ 
making Ji 2 non-symmetric under the exchange of the 
two qubits). It is striking that we also find in our model 
that C{t) — Cp{t} for specific initial states even though 
the dynamics in the absence of measurements - and thus 
the density matrix concurrence Cp(() - arc not the same 
in the two models (here the two qubits are coupled to a 
common bath, whereas they are coupled to distinct baths 
in Ref. (H). 

VIII. CONCLUSION 

We have found explicit formulas for the mean concur- 
rence C{t) of quantum trajectories and have shown that 
the measurements on the baths may be used to protect 
the entanglement of two qubits. These results shed new 
light on the phenomenon of entanglement sudden death. 
For independent baths, C{t) is either constant in time 
or vanishes exponentially with a rate depending on the 
measurement scheme only, whereas for a common bath 



C{t) depends strongly on the initial state and may co- 
incide with the concurrence Cp(() of the density matrix 
for some initial states. A constant C{t) implies a perfect 
protection of maximally entangled states for all quan- 
tum trajectories. In the case of pure dephasing and for 
Jaynes-Cumming couplings at infinite temperature, we 
have found measurement schemes independent of the ini- 
tial state of the qubits which lead to such a perfect en- 
tanglement protection. Despite obvious analogies, this 
way to protect entanglement differs from the stra tegy 
based on the quantum Zcno effect proposed in Ref. |l5| . 
In fact, in the QJ and QSD models considered here the 
time interval between consecutive measurements is not 
arbitrarily small with respect to the damping constant 
7^^. In the QJ model this time interval dt must be cho- 
sen such that the jump probability dp{t) cx 7 dt is very 
small but one cannot let 7 dt go to zero since this would 
amount to replacing dp{t) by and e"'^"'"'"^* by e"'-'^"'^* 
in ([6]). In contrast, a perfect entanglement protection is 
reached in [l5j in the idealized limit 7dt ^ (i.e., when 
the measurements completely prevent the decay of the 
superradiant state HI]). 

For independent baths, C{t) is strictly greater than 
Cp(j) if the latter concurrence vanishes after a finite time. 
Therefore, if there exists a measurement scheme such 
that the mean entanglement of formation E(t) is equal 
to the entanglement of formation of the density matrix 
(which would imply C'(t) < Cp(^t))i this scheme must nec- 
essarily involve measurements of non-local (joint) observ- 
ables of the two baths. Let us finally note that it should 
be possible to check our findings experimentally by using 
similar optical devices as in Ref. [J]. 
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